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We derive from a microscopic Hamiltonian a set of stochastic equations of motion for a system 
of spinless charged particles in an electromagnetic (EM) field based on a consistent application of 
a dimensionful 1/c expansion of quantum electrodynamics (QED). All relativistic corrections up 
to order 1/c 3 are captured by the dynamics, which includes electrostatic interactions (Coulomb), 
magnetostatic backreaction (Biot-Savart), dissipative backreaction (Abraham-Lorentz) and quan- 
tum field fluctuations at zero and finite temperatures. With self-consistent backreaction of the EM 
field included we show that this approach yields causal and runaway-free equations of motion, pro- 
vides new insights into charged particle backreaction, and naturally leads to equations consistent 
with the (classical) Darwin Hamiltonian and has quantum operator ordering consistent with the 
Breit Hamiltonian. To order 1/c 3 the approach leads to a nonstandard mass renormalization which 
is associated with magnetostatic self-interactions, and no cutoff is required to prevent runaways. 
Our new results also show that the pathologies of the standard Abraham-Lorentz equations can be 
seen as a consequence of applying an inconsistent (i.e. incomplete, mixed-order) expansion in 1/c, 
if, from the start, the analysis is viewed as generating a low-energy effective theory rather than an 
exact solution. Finally, we show that the 1/c expansion within a Hamiltonian framework yields 
well-behaved noise and dissipation, in addition to the multiple-particle interactions. 



I. INTRODUCTION 

The backreaction of a charged particle interacting 
with the electromagnetic field involves a number of fa- 
mous problems including acausality, in the form of pre- 
acceleration, and runaway solutions to the Abraham- 
Lorentz equation [ ]. A number of approaches to re- 
solving these problems have been developed, including 
replacing point particles with extended objects [2-4], 
treating the electromagnetic field interaction perturba- 
tively and truncating at a specific order [5], or replacing 
time-local differential equations of motion with nonlocal 
integro-differential equations of motion [G, 7]. Yet many 
problems related to the backreaction remain poorly un- 
derstood or unresolved. The interaction of a local (par- 
ticle) degree of freedom with a nonlocal environment (a 
field) can also result in surprising dynamics, for example 
the vanishing of the radiation-reaction force on a uni- 
formly accelerating charge. 
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To find the backreaction on a single charged particle 
within an open systems framework, one integrates out 
the field degrees of freedom following well-known proce- 
dures [6]. Starting from a nonrelativistic Hamiltonian of 
a charged particle coupled to an electromagnetic field, the 
procedure gives the Abraham-Lorentz Langevin equation 
for a structureless point charge coupled to the electro- 
magnetic field [5, 6, 8]. fncluding driving by external 
forces, the stochastic equations of motion are 

m x(t) = F oxt (t) + 2e 2 7 Qx(t) - e|(f) , (1) 

backreaction noise 

where 70 = l/12-7T£oc 3 and l-(t) is quantum-field induced 
noise. These equations of motion can describe either 
classical trajectories or Heisenberg-picture operators. In 
the latter case, the induced noise is operator valued and 
has both a complex-valued noise correlation and a state- 
independent commutator. As noted above, Eq. (1) is not 
manifestly causal and exhibits runaway solutions. The 
result in Eq. (1) is essentially equivalent to backreaction 
with a supra-Ohmic bath in quantum Brownian motion 
(QBM) [5, 9]. However, relativistic QED has nonlinear 
coupling of system and environment, whereas QBM has 
bilinear coupling. Consequently, the physics of QED is 
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far more subtle and complex. 

In this paper, we construct an alternative derivation of 
the stochastic equations of motion for a system of spin- 
less charged particles based upon a consistent application 
of a dimensionful 1/c expansion. We show that this ap- 
proach yields causal (i.e., no-preacceleration) and stable 
(i.e., runaway free) backreaction, provides new insights 
into previous results, and naturally leads to equations 
with a form consistent with the (classical) Darwin Hamil- 
tonian and operator ordering consistent with the Breit 
Hamiltonian. The role of magnetostatic interactions in 
both mass renormalization and as a form of dissipation- 
less backreaction is also revealed. 1 For a single particle, 
our damping admits correspondence with the casual and 
stable equations of motion obtained previously by Ford 
and O'Connell [ ]. 

To order 1/c 3 the approach leads to a nonstandard 
mass renormalization which is associated with magneto- 
static self-interactions, and (to this order) no cutoff is 
required to prevent runaways. An important question is 
whether the 1/c expansion automatically generates con- 
sistent low-energy behavior that is always insensitive to 
particle structure. Although it will be complicated, our 
method provides a systematic method to consistently ex- 
tend the analysis to higher orders for a closer investi- 
gation of this question. This is important because the 
standard analysis giving the Abraham-Lorentz equation, 
viewed perturbatively in 1/c, is of mixed order. Our re- 
sults show that inclusion of incomplete higher-order in- 
formation (that is, some but not all terms beyond 1/c 3 ) 
can be identified as the cause of the Abraham-Lorentz 
equation's perturbative instability, from an effective the- 
ory perspective. Equivalently, our results indicate that 
in terms of a 1/c expansion the dipole approximation is 
inconsistently applied in the standard derivations, and 
that if all terms present in the standard calculation are 
to be included, then additional multipolc field corrections 
will be also required to the same order. 

Our analysis of the single particle dynamics also high- 
lights a subtle but important connection between dissi- 
pation and noise. There are two standard calculations 
for radiation reaction which arise from different choices 
of gauge: Coulomb or electric dipole. We show that the 
former (see, e.g., [8] for an example) has a problematic 
description of noise: it is not manifestly thermal in the 
sense that the stochastic process present in the Langevin 
equation is not sampled from an independent thermal dis- 
tribution constructed from the field's Hamiltonian acting 
as a reservoir. The latter (see, e.g., [5] for an exam- 
ple) has ordinary noise but involves integration kernels 
which are approximately given by the second derivatives 



of a delta function, making them more pathological than 
the those in the Coulomb gauge. In other words, either 
the noise or dissipation has undesirable features in the 
usual derivations. We show that the 1/c expansion in 
the Coulomb gauge, with a correct delineation of system 
and environment, naturally avoids these problems, yield- 
ing well-behaved noise and dissipation in addition to the 
multiparticle interactions within a unified framework. 

In Sec. II, we first analyze nonrelativistic quantum par- 
ticles in the electromagnetic field starting from the usual 
Hamiltonian but instead of taking the dipole approxi- 
mation we expand in orders of powers in 1/c. Contin- 
uing with the usual calculation we reproduce the stan- 
dard results, but eventually see that they are "mixed 
order" in powers of 1/c. In Sec. IV, we work from the 
Coulomb gauge and throughout the calculation we con- 
sistently remove all terms in the open-system dynamics 
beyond C(l/c 3 ) . The usual formulations can in principle 
be applied but it is not as transparent and is more com- 
plicated than necessary, as unneeded higher order terms 
are first inconsistently kept and only later removed. Sec- 
tion IV also shows that the order 1/c 3 equations of mo- 
tion yield stable backreaction and self-consistent noise. 
Section V discusses the regime of validity of our analysis 
and discusses future directions. In App. A, we review 
the derivation for backreaction in QBM which parallels 
our analysis in many respects. Other characteristics of 
the QBM model which are used here, such as the renor- 
malization and integration kernels, are more thoroughly 
discussed in [ ] . In Appendix B we solve the field equa- 
tions of motion sourced by the particles, without invoking 
the dipole approximation. 



II. STANDARD "NONRELATIVISTIC 
RADIATION REACTION 



We first review the derivation of the standard Hamilto- 
nian for the motion of a system of N particles, while also 
defining the notation we will use in the rest of the paper. 
We start with the Lagrangian £ sys for the motion of the 
system defined as the collection of particles j = 1,...,N 
with mass rrij and charge ej for the j th particle and the 
Lagrangian £ cnv for an electro (E)- magnetic (B) (EM) 
field acting as its environment: 



By electrostatic and magnetostatic we refer to the lowest-order 
(in 1/c) electric and magnetic fields. These fields are instanta- 
neous and accompany even static charge and current distribu- 
tions, though our system is not static. The higher-order rela- 
tivistic contributions contain retarded and radiative effects. 



•^QED — £sys + Ant + -Ccnv , (2) 

£ sys = \ m J ±2 j ' ( 3 ) 
3 

C cm = I dV^ {E(x') 2 - c 2 B(x') 2 } • (4) 
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The interaction between the system of charged particles 
and the EM field is given by 

Ant = J d 3 x' {J(x, x') • A(x') - p(x, x') «£(x')} , (5) 
= 5^e J -{x r A(x i )-0(x J -)} , (6) 



where p and J are the charge and current density cou- 
pled to the scalar and vector potentials <fi, A respectively 
which are related to the electromagnetic fields by 



E(x) = - V0(x) - ^A(x) , 
B(x) = VxA(x). 



(7) 
(8) 



Expressing the vector potential in terms of the spatial 
Fourier modes with wavevector k and polarization e k 
gives 



AW= (2^72 /^ 3 fcE A M,(x) 



(27T) 

A k , £fc (x) = { e +* kx a k , £fc + e- k x a* 



(9) 



4' (io) 



where u>k = ck. To satisfy the commutation relations, 
the conjugate momentum of the field is then given by 



7r W = /^£7T k , efc (x) 



(11) 



7Tk, efc (x) 



^2/e u> k 



zk-x 



at. 



tkx t 



(12) 



To systematically treat the position dependence of the 
coupling, we write 

Ak,e fc (x) = cos(k • x) A k efc - sin(k ■ x) , (13) 

7T k , £(c (x) = cos(k • x) 7T k , £(j + sin(k • x) e uj k A k , £fc , 

(14) 



where 



A k , £fc = A k . £fc (0), 
7Tk, £fc = 7T k , £fc (0), 



(15) 
(16) 



are x-independent field operators that more closely cor- 
respond to the "positions" and "momenta", g k and 7r k , 
of the reservoir oscillators for QBM, defined in App. A. 



A. Coulomb-Gauge Hamiltonian 

Choosing the Coulomb gauge V • A = 0, the scalar 
potential 4> becomes the nondynamical Coulomb poten- 
tial and the vector potential A is purely transverse, i.e. 



k ■ e fc = 0. The system plus environment Hamiltonian 
may then be expressed as 



h qed - E 



P './ A ( x ./r 

2m a 



(17) 



Hfidd — Id 3 



J d 3 x' 1 - jV7T(x') 2 + e Q c 2 [VxA(x')f) 



(18) 



Expanding in field modes, the Hamiltonian is 

H ncld EE Jd 3 kJ2\ {eo 1 nl,e h + ^le k } • (19) 

The Heisenberg equations of motion for the system as 
driven by the field are then 



TO, TO, 



-A(xj-) , 



(20) 



P = -ejVj^x,-) + —Vj [ Pj - < , A(x,)] 2 , (21) 

and the Heisenberg equations of motion for the environ- 
ment as driven by the system are given by 
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— ^IWk.x,),^}, (22) 



A ki£fe — "TTk^fc . 

J 

7Tk, e , = -eo^fc A k , £fc + y {cos(k • Xj), ±j} , (23) 

where {A, B} = A B + B A is the anticommutator. From 
App. B, the driven solution can be expressed in the 
manifestly-Hermitian form 

A( Xi ,t) = (*)-5> {(m£ ***)(*) + (/*« *ii) f (*)} . 

( 24 ) 

with the convolutions defined as 

(/4 /" dt' / x^[x i (t),x i (t / );t-* , ]x i (t / ), 

Jo 

where the dissipation kernel fi A is defined exactly in 
App. B and resolved (to lowest orders in 1/c) in Sec. III. 
The dissipation-kernel integral arises from the system 
driving the field, whereas the operator ^ i (t) corresponds 
to the homogeneous evolution of the field. 

Multiparticlc Langevin equations of motion can then 
be obtained by substituting Eq. (24) for A back into 
the system equations of motion (20)-(21). We will not 
write down the multiparticle equations in this approach, 
however, since they are more cumbersome than the form 
we obtain in Sec. IV. We analyze the single-particle case 
in Sec. IIIC. This gives the Abraham-Lorentz Langevin 
result in [ ]. 

Although the Langevin equation derived from Eq. (24) 
is valid as an equation of motion for the Heisenberg oper- 
ators, proper global boundary conditions must be applied 
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so that the homogeneous-evolution operator %f(t) can be 
given the interpretation of independently-sampled ther- 
mal noise. (For the simpler case of QBM these details 
are described in App. A3-A4.) To obtain statistically- 
independent thermal noise, the Langevin equation must 
be derived assuming an initially factorized system and 
environment state p sys+cnv = p sys ® p onv with the en- 
vironment initially in its equilibrium state. (Of course 
Psys+cnv ls n °t an equilibrium state at the initial time, 
only p cnv .) However the field equations in Eq. (24) are 
expressed in terms of the particle velocity Xj(t) rather 
than the canonical momentum pj (t) and this implies that 
the initial state of the system plus environment cannot 
immediately be placed into the standard form with re- 
spect to the Langevin equation. 

To see this, note that for a Langevin equation sourced 
with x,(t) we must supply initial data Xj(0), Xj(0), etc., 
for the particle that is independent of (uncorrelated with) 
the initial state of the environment. Most naively, such 
an initial state would have the factorized form 

Psys+cnv = Psys( X J x ) ® Pcnv(A, T) ■ 

However, because the velocity Xj(t) does not commute 
with the canonical momentum 7T, such a state with 
p env (A, 7T) cannot represent a (field) equilibrium state 
with respect to Hamiltonian (19). It is also likely that 
such a "factorized" initial state is not even a proper den- 
sity matrix as the product of two noncommuting positive- 
definite matrices is not necessarily positive definite. A 
nontrivial initial state that did take the above form 
would necessarily involve, therefore, an initial nonequi- 
librium state of the environment that implicitly depends 
on p sys (x,x). The operator- valued quantum noise £j (t) 
would then be sampled from this initial nonequilibrium 
state of the environment and consequently it would not 
represent standard thermal noise. This would also nec- 
essarily imply that the initial state of the system would 
always be dependent upon the realization of the noise. 
To obtain standard (initially uncorrelated) noise from an 
initially equilibrium environment, the initial state of the 
system must be of the form p sys (x, p), represented in 
terms of canonical coordinates x and p. 



B. Electric-Dipole Gauge Hamiltonian 

To obtain a Langevin equation with standard, statisti- 
cally independent noise, we can make the canonical trans- 
formation 



ex • A(x) = ex • A(x) — e ^[ x ■ A(x)] , 



and neglect the total derivative in the action (see, e.g. 
[ ]). The Hamiltonian becomes 

H QED = E|^-+ e ^(^)}+ ( 25 ) 

/ dV { 2^ [7r(x ' } ~ P(X ' X ' )X]2 + [V X A(X ' )]2 } ' 

which is analogous to the QBM Hamiltonian (A6). The 
interaction is 

Hint = ~y] ej xj ■ 7T(Xj-) , 

3 

after the p 2 x 2 term is included in the renormalization of 
the potential <j>. 

The equations of motion for the system driven by the 
field are 

x,.= Pi, (26) 
rrij 

Pj = -eiV^b a re(Xj) + ejVj [xj ■ 7T(xj)} , (27) 

where we take the bare potential to include the diver- 
gent contribution from p 2 x 2 in the Hamiltonian. The 
analogous analysis for QBM is described in App. A. 

The equations of motion for the field driven by the 
particle arc 

A k , £fc = — ^M* - e 3 cos(k • x 3 -) Xj , (28) 

3 

7Tk, Efc = -e ul A kj£fc + e u> k ^ e 3 sin(k • x,-) Xj . (29) 

3 

As calculated in App. B, the driven solution expressed in 
the manifestly-Hermitian form is 

TT(x^) = £f {(/ig *x i )(*) + (M5 *x^)t(t)} , 

(30) 

where the dissipation kernel fifj is resolved (to lowest 
orders in 1/c) in Sec. III. 

In parallel to the analysis of the previous section, mul- 
tiparticle Langevin equations of motion can be obtained 
by substituting Eq. (30) for 7T into the system equations 
of motion (26)-(27). Again, we delay writing them down 
because in this approach the multiple-particle equations 
of motion are more cumbersome that the form given in 
Section IV. The single-particle case, however, is analyzed 
in Sec. IIIC. 

Unlike the analysis from the previous subsection, the 
field variables are now expressed in terms of the canon- 
ical variable Xj, in contrast to Xj, and consequently 
the stochastic variable £f (i) can be interpreted as 
statistically-independent noise . Issues due to noncom- 



2 In fact, £f(t) can be (perturbatively) identified with a stochastic 
electric field, though 7T cannot be identified with the electric 
field as it does not contain the electrostatic fields and necessarily 
contains the magnctostatic fields. 
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mutativity of system (particle) and environment (field) 
coordinates discussed in Sec. II A do not occur in this 
gauge since an initial state of the form 



P BY s( x >P)®Penv( A > 7r ) 



(31) 



is consistent with the required initial data for particle po- 
sition in the Langevin equation, with p env (A,7T) at the 
same time an equilibrium state of the field Hamiltonian. 
As will be reviewed in the next section, however, the 
damping kernel for fifj is more pathological than for fifj . 
The latter is approximately a delta function (Ohmic), 
whereas the former is approximately the second deriva- 
tive of a delta function (supra-Ohmic). 

To summarize, stochastic equations for motion ob- 
tained for the Coulomb gauge with x ■ A coupling in the 
particle-field interaction have noise that is not statisti- 
cally independent of the particle's requisite initial data, 
and this will lead to severe complications in both the 
interpretation and evaluation of the resulting Langevin 
equations. In contrast, stochastic equations of motion for 
the electric-dipole gauge with x-7T coupling have statisti- 
cally independent noise, but more pathological damping, 
as described in the next section. 



III. ELECTROMAGNETIC DAMPING 
KERNELS 

We next analyze the nonlocal integration kernels which 
arise for the Langevin equations. It is useful to first define 
the commutators of vector operators 

[X,Y] = XY T -YX T . 

This object is a matrix whose entries are ordinary com- 
mutators. In this notation, the dissipation kernel associ- 
ated with A-coupling from Eq. (24) is exactly specified 
in App. B and is approximately given by the field com- 
mutator 



Az A [x(t),y(*');M'] 



1 

2ih 



[A(x(i),i),A(y(0,i') 



with phase discrepancies of the order O (fi 2 /c 3 ). These 
phase discrepancies denote quantum-relativistic devia- 
tions from the quasistatic approximation wherein we take 
x(i) and y(t') to be adiabatic variables inside the field 
correlations. This is largely equivalent to the dipole ap- 
proximation for one particle, but it allows us to system- 
ically extend the analysis to the multiparticle case while 
keeping tracking of the order of the approximation in 
powers of 1/c. More details on this relativistic "multi- 
pole" expansion are given in App. B. 

Let us denote the time dependence in the field opera- 
tors which arises from x(i) as intrinsic and the explicit 
time dependence in A(r,t) at some fixed location r as 
extrinsic. With this distinction, the dissipation kernel is 
stationary with regard to extrinsic time dependence, i.e., 

MaW*). y(*0; *> = maK*), y(t');t-t'] . 



It is not stationary, however, with regard to the intrinsic 
time dependence of x(t) and y(t'). The quantum dissipa- 
tion kernel, however, is approximately spatially station- 
ary in the sense that 

H A [x(t),y(t');t,t'} « n A [ X (t)-y(t');t-t'} , 

with O (h 2 /c 3 ) phase discrepancies as discussed in 
App. B. 

The positive-definite and Hermitian damping kernel is 
then given by 



d 



(32) 



/x[x(t),y(i');t,i'] = -^^[xW,y(t');M'] , 
A[x(t),y(0;w] = iWy[x.(t),y(f)M . ( 33 ) 

where the partial derivative and Fourier transform ne- 
glect any intrinsic time dependence in x(t) and y(t'), 
similar to the analysis in [11]. As discussed in App. B, 
we may neglect the time dependence intrinsic to y(t') 
when taking a total derivative, as the corrections are of 
higher order in v/c. 

Fourier transforming with respect to the extrinsic time 
variables, the classical or quasistatic A-coupling damping 
kernel (without cutoff) is 



7o 



27o{Si( 



rw 




roj 






)+S ( 




) 


c 


c 





1 



12ire c 3 ' 
in terms of the functions 
3(z 2 - 



1) sin(z) + z cos(z) 



Si(z) 

- 3 (z 2 — 3) sin(z) + 3 z cos(z) 

°0(Z) = - o -5 



(34) 
(35) 

(36) 
(37) 



A cutoff regulator \ can be easily inserted by multiplying 
the right-hand-side of Eq. (34) with a function 



X (lo/A) : [0,co)^ [1,0), 



(38) 



that vanishes sufficiently fast. 

In the coincidence limit we recover the usual Ohmic 
damping: 



lim n 7A[ r ^} = 2 7o , 
lim7i[r;t] = 2-y <5(*) 



(39) 
(40) 



and for particles that are far separated all cross correla- 
tions vanish: 



lim 7 A [r;cj] =0, 

r— >oo 

lim 7a [r; t] = . 



(41) 
(42) 



In Fig. 1 we compare these special functions to sinc(z), 
which is the result obtained from the simpler but anal- 
ogous case of coupling to a scalar field. In Fig. 2 we 
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FIG. 1. Comparison of sine (bold), Si, and So (dashed). Sine 
and Si are extremely qualitatively similar, both being unity 
at zero whereas So vanishes at zero. 
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FIG. 2. The same functions as in Fig. 1, but in the time do- 
main: the rectilinear distribution (bold), So, and Si (dashed). 

compare these functions in the time domain. The fre- 
quency domain is useful for noting the Markovian limit at 
no = and the decorrelation which occurs for rw — > oo. 
The time domain is useful for illustrating that the damp- 
ing kernel is (perturbatively) causal and constrained to 
the light cone. As is well known, this latter property is 
not true of the quantum noise kernel. 

A. Dissipation and Damping with A-coupling 

To describe quasirelativistic damping, we can inte- 
grate by parts the A-coupling dissipation kernel to obtain 
terms corresponding to damping, mass renormalization, 
and dissipationless backreaction. In the quasirelativistic 
regime we treat the positions x as quasistatic within the 
convolution (as discussed in App. B) 

Q4*z)(t)= f dt'ix A [xi(t)-xj(t');t-t']z(t'), 
Jo 

where z(t) is some arbitrary source. Integrating by parts, 
and neglecting the intrinsic time dependence in the po- 



sitions (which is equivalent to ignoring higher-order cor- 
rections in v/c), gives 

(My **)(*) - (7«**)(*) - TA[xi(*)-Xi(t);0] *(t) 
+ 7 A [x i (t)-x*(0);i]z(0) > 

where the final term is a transient slip. Although the slip 
evolution arises during the same period of time as the 
conventional Abraham-Lorentz acausality, its role is well 
understood from QBM studies where it has been shown 
to be unproblematic (see App. A and Ref. [10] where 
the initial evolution for initial factorized system plus en- 
vironment states in terms of so-called "slip" and "jolt" 
evolution is reviewed). It may therefore be discarded for 
our purposes; doing so is equivalent to choosing an initial 
state which is more properly correlated to the environ- 
ment. The difference between QBM and QED is that 
with the relativistic field, not only are there instanta- 
neous self-transients from J^t), there are also retarded 
cross-transients from 7^(t). Figure 2 shows that the 
electrodynamic slip only produces a transient effect for 
t < r/c. Essentially, for factorized states of the particles 
and field, the particles are completely unaware of each 
other's existence until photons travel between them and 
suddenly (and violently) establish correlations at their 
first mediation. 

Unlike QBM, ~f A [rij(t); 0] contains both self interac- 
tions, which renormalize the mass, and field-mediated 
interactions. Using Eq. (34), this term becomes 

3 l + ff T 
7aM] = -c l0 |r| , 

and we see the emergence of the Darwin Hamiltonian 
(e.g. see Jackson [1] Sec. 12.7). The magnetostatic en- 
ergy of the Darwin Hamiltonian (which is classical) is 
given by 

v, ; p,'7 ( !•„.:«) '; p.. 

3 3 

The quantum extension of the Darwin Hamiltonian re- 
quires an operator-ordering prescription that our analysis 
will provide upon applying these results to the Langevin 
equation. 

These results exemplify how the equations of motion 
that result from an A-coupling interaction give consistent 
formulas for damping and backreaction, despite the noise 
being problematic as previously discussed. 



B. Dissipation and Damping with Electric-Dipole 
Coupling 

In this subsection we describe the sense in which the 
electric-dipole interaction gives a less desirable descrip- 
tion of damping and backreaction, despite providing a 
consistent description of thermal noise, as previously de- 
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scribed. For 7T-coupling to the field, the dissipation ker- 
nel is approximately given by 



1 

2th 



[7r(x(t),t),7r(y(t'),t' 



with the same O (h 2 /c 3 ) phase discrepancies discussed in 
App. B. The 7T-coupling dissipation kernel is related to 
the A-coupling dissipation kernel via 



d 2 

7 7r [r;t,i] =__ T4 [ r> t,t] ) 



(43) 
(44) 



This result is equivalent to dissipation from a supra- 
Ohmic bath in QBM. As an integration kernel 7^ is 
relatively pathological and must be integrated by parts 
twice to obtain the well-behaved kernel 7^4. Such an in- 
tegration is straightforward for the single-particle theory 
in the dipole approximation, but for the multiparticle 
and higher-order relativistic theory this generates many 
additional terms and limits which must be carefully ana- 
lyzed in constructing the Langevin equation. Although, 
with care, it should be possible to proceed this way, the 
approach in the next section is clearer and more straight- 
forward. 



C. Standard Dipole Calculation 

Before presenting our effective equations of motion in 
Sec. IV, we first review the derivation of the standard 
results for a single particle. Taking the dipole approxi- 
mation we drop all position dependence in the integration 
kernels. Substituting Eq. (30) (integrated by parts once) 
into Eqs. (26)-(27) we obtain the Langevin equation 

mx(f ) = -eV0(x) - 2 e 2 (7* * x) (t) +e£ w (t) , (45) 

discarding the transient terms. This equation of mo- 
tion is runaway free for bare mass m > as reviewed 
in App. A 5. 

Integrating-by-parts the 7r-damping integral in the 
Langevin equation (45) two times, to obtain A-coupling 
damping, discarding additional transient terms, and 
grouping like terms, we obtain the standard Abraham- 
Lorentz-Langevin equation 

m ren x(t) = -eV0(x) + 2e 2 ( 7 A * x)(i) - e £ A (t) , (46) 



= m + 2e 2 7 A (0), 



(47) 



where m rcn is the renormalized particle mass. The same 
result is obtained by integrating-by-parts the dissipa- 
tion integral in the Langevin equations (20)-(21). Note 
that (7 A * x)(i) rj Yo'x(£) in the high cutoff ("point- 
particle") limit. Positive bare mass to, which is re- 
quired for runaway-free motion, requires a finite cutoff 
A in the field modes such that 2e 2 j A (0) < m ren , not- 
ing that 7 (0) oc 70 A. Therefore positive bare mass 



requires an ultraviolet cutoff A, or equivalently a form 
factor that cuts off the particle field coupling, on the or- 
der of £oC 3 TO ren /e 2 or smaller. This is directly seen from 
our stability analysis in App. A 5, or simply by noting 
that with a negative bare mass the system + environ- 
ment Hamiltonian no longer has a lower bound in its 
energy spectrum. 



IV. 1/c EXPANSION OF QED 

In the previous section, we contrasted how the stan- 
dard results are obtained when the interaction is ex- 
pressed in terms of either x • A or x ■ 7T coupling. We 
have highlighted the problem of simultaneously obtain- 
ing well-behaved noise and backreaction. In [7], Ford 
and O'Connell showed that the equations of motion for a 
charged particle with structure can be order-reduced to 
obtain the structure-independent equation of motion 



1 



2e 2 7o d 
to dt 



F ext (x) - e£(t) 



which is accurate to lowest order in the timescale 



e 2 /67T£oTOc 3 ~ ah/mc 2 



(48) 



(49) 



where a — e 2 /A-keqHc is the fine-structure constant. For 
bound states, this is equivalent to the dimensionless order 
aAE/mc 2 , where AE denotes the relevant energy-level 
transitions of the system. For driving forces, the expan- 
sion parameter is a hu /mc 2 , given a driving frequency 
uj. In either case the result is of order 1/c 3 , and the e 3 £ 
term is negligible given the approximations which have 
been made in constructing this equation. 

Eq. (48) is causal, runaway-free, and is in general dis- 
sipative (i.e. accelerated motion is damped). The Ford 
and O'Connell result is generic in the sense that any 
particle, with reasonable assumptions about its struc- 
ture (form factor), will have backreaction of this form 
in the weak-backreaction limit. In this sense, the result 
is the universal (effective) low-energy, and equivalently 
lowest order in 1/c, result for backreaction. Below, we 
obtain stochastic equations of motion with damping con- 
sistent with Eq. (48), but from a different approach that 
provides additional insights into the nature of backreac- 
tion. Whereas the more-standard approach of Ford and 
O'Connell considers Eq. (1) to be of order C(l/c 3 ) and 
their Eq. (48) to be order reduced, our approach will 
generate equivalent damping to Eq. (48) which is strictly 
C(l/c 3 ), and not order reduced. Within our analysis, 
the standard Eq. (1) is mixed-ordered, and this is the 
primary source of its pathological behaviors. An advan- 
tage of our method is that it allows for a consistent and 
systematic extension to higher orders which includes both 
higher-order relativistic corrections, and higher-order ef- 
fects of backreaction. 

Our approach fundamentally relies upon consistent ap- 
plication of a dimensionful 1/c expansion around the 
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Hamiltonian 



H 



NR — 



t-r 1 2m, 



4-?. 47re |r< 

i<3 



which is viewed as the lowest-order, "nonrelativistic" 
Hamiltonian in this approach. All of the additional terms 
that arise in the QED Hamiltonian are then viewed as 
perturbations, whose order is based on the powers of 1/c. 

We begin with the A-coupling, Coulomb-gauge Hamil- 
tonian in (17). From the previous analysis we observe 
that the self-field generates, at lowest order C(l/c 2 ), 
magnetostatic and renormalization terms (see Sec. Ill A). 
Because our approximation is based on a systematic ex- 
pansion in powers of 1/c, and here we work to 0(l/c 3 ), 
we effectively drop the A 2 system-field interaction terms 
in the Hamiltonian.' 5 These neglected terms are all of 
order (v/c) 4 , a (AE /mc 2 ) (v / c) 2 , and so, making them 
C(l/c ). For now we assume the usual nonrelativistic ki- 
netic energy for the particles, but for consistency we in- 
clude the relativistic corrections at the appropriate order 
in Sec. IV C. Therefore, the consistent 0(l/c 3 ) effective 
Hamiltonian in the Coulomb gauge is 



(50) 



V { Pl ■ A(X() + A(x<) ■ p 4 } + H ficlcl + G(l/c 4 ) 



The Heisenberg equations of motion for the system 
driven by the field for the effective Hamiltonian Heqed 
are 



x.j — A(x.j) , 

m, rrij 



(51) 



p., = -eiVi(/>(xi) + w^-Vi {pi ■ A(xi) + A(xi) • pj , 
zm, 

(52) 

and the Heisenberg equations of motion for the environ- 
ment driven by the system are 

A k , efc = ^7T k , £fc +- ] —Y J ^t {sin(k ' ' 

e e uj k ^ Zm., 

3 J 

(53) 

7Tk,e k = -e Q u 2 k A k , efc + {c° s ( k • x j), Pj} ■ 



(54) 



As derived in App. B, the driven solution expressed in 
manifestly-Hcrmitian form is 



A( Xi ,t) = ± {04 * p,o(*) + 04 * p,) f w} 

3 3 

(55) 



3 For external fields this approximation should not be made, and 
any A cx t should be included by appropriately translating the 
system momenta. 



where the £,(i) are uncorrelated processes for i ^ j. Fol- 
lowing the approach in Sec. Ill A, integrating the dissi- 
pation kernel by parts gives 



A[xi(t),t] =£i(*)-2— 7oPi (t) 



(56) 



J2 { -r[Tij(t);0]—PiV) + (7[r«(*);0]-^Pj(t) 

'— J TOo V TOo 



Here we have neglected the transient slip terms (which, 
as noted before and as is reviewed in the appendix, can 
be verified to not induce runaways or acausal behavior) 
and corrections of higher order in v j c. 

Substituting the driven field solution into the sys- 
tem equations of motion gives the renormalized multiple- 
particle Langevin equation 



Pi gj 

TOo m. 



x , = ; £ i + 2 -% 70 Pi + V Pl Y *3 , ( 57 ) 



pi = - eiV Xl (f>(xi) - V Xi ^ Vij , 

3¥=i 



(58) 



Details of the renormalization will be explained further 
in Sec. IV B. The quantum magnetostatic potential V is 
given by 

V « = [Pi"T[r«;0] Pi + P.S r.„:0 Pi ] 

~ \ S ^ TTx [p * P ^ T l[VlJ ; 0] + l[ri3 ; 0] Pl p ? ] ' 

(59) 

The last two terms have been expressed with a spatial 
trace 



Tr a [M] = M x 



M, 



yy 



M zz 



(60) 



to keep the proper Hilbert-space operator ordering. 
The complexity of this expression is due to the non- 
commutativity from both Hilbert-space and spatial (3- 
vector and 3x3 matrix) operations. This quantum mag- 
netostatic potential is consistent with the (classical) Dar- 
win magnetostatic potential (also see Sec. Ill A) 



V ij =-2^p? 7 [r ii ;0; 



r nl 3 1 
7[r;0] = ^70- 



r r 



Pj 



(61) 
(62) 



Moreover, the operator ordering in Eq. (59) is consistent 
with the Breit equation for spin-1/2 particles [12] (often 
named "orbit-orbit interaction" in that context), when 
one ignores all of the Pauli spin matrices. 

Equations (57)-(58) are one of the main results of this 
paper. They are multiparticle stochastic equations of mo- 
tion consistent through order C(l/c 3 ) in the field influ- 
ences with, as we will show, stable (runaway-free) back- 
reaction even as the field cutoff goes to infinity. [O (l/c 3 ) 
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kinematics will be given in Sec. IV C] Notice that our re- 
sults show that the mass renormalization associated with 
the dissipative backreaction, at C(l/c 3 ) , is due to mag- 
netostatic self-interactions. Our effective Hamiltonian 
treatment also correctly produces the second-order mag- 
netostatic corrections without extraneous fourth-order 
terms (as Breit accidentally included in his first calcula- 
tion [12, 13]). This result is not present in the standard 
treatments. 

Our analysis suggests that the pathologies of the stan- 
dard Abraham-Lorentz equations can be viewed, from 
the perspective of a 1/c effective theory expansion, 
as a consequence of performing a mixed-order calcula- 
tion. In terms of a 1/c expansion starting from the 
Coulomb-gauge Hamiltonian, there are C(l/c 4 ) multi- 
poZe-expansion terms in the p • A(x) interaction terms 
that are of the same order as the lowest-order dipole 
terms present in the A(x) 2 interaction. This reveals that 
uniform application of the dipole approximation by it- 
self to the Coulomb gauge Hamiltonian is inconsistent 
with the 1/c expansion, in the sense that it keeps some 
terms of C(l/c 4 ), but discards others. For a consistent 
C(l/c 3 ) expansion, effectively the dipole approximation 
to the p • A(x) interaction is kept and the A(x) 2 inter- 
action is dropped entirely, as we have done. 

Historically, the Abraham-Lorentz equation has most 
often not been viewed perturbatively, but instead as an 
effort to obtain an exact, nonrelativistic point-particle 
result. Our philosophy has been to instead look for a 
consistent low-energy, effective theory from the begin- 
ning. Giving up any claim to an exact solution, we gain 
new insight into how a perturbatively consistent solution 
arises within a framework that can be systematically ex- 
tended to higher orders. An alternative approach is to ap- 
ply order reduction to the standard result, which can be 
iterated into the Ford-O'Connell equation (48), with ad- 
ditional backreaction contributions of higher order in 1/c 
(and higher-order derivatives). One then discards these 
higher-order terms in the name of order reduction. Com- 
bining our relations (57) and (58) reproduces the same 
damping as the Ford-O'Connell equation, however, with- 
out any higher-order terms to be order reduced. Thus, 
in terms of a 1/c expansion, the damping in the Ford- 
O'Connell equation is strictly 1/c 3 (and thus consistent), 
whereas the pathological backreaction in the Abraham- 
Lorentz equation is of mixed order. There is an interest- 
ing parallel here to Breit's original calculation [13], where 
pathological equations resulted when Breit accidentally 
included a fourth-order operator in a calculation which 
was only accurate to second order. 

A. Noise and backreaction stability 

As previously discussed in Section II A the "noise" 
for the standard v • A coupling Langevin equation can- 
not be independently sampled from a thermal distribu- 
tion. The standard calculation only gives well-behaved 



noise within the electric-dipole gauge. In contrast, our 
analysis, to order C(l/c 3 ) starting from the Coulomb 
gauge, gives well-behaved noise and dissipation. Compar- 
ing the derivations in Sections II A and IV, we see that 
our noise and the standard electric-dipole gauge noise 
only differ by contributions from the A(x) 2 interaction 
of the order 1/c and higher, or more specifically (v/c) , 
a (AE/mc 2 ) (v/c) 2 , etc.. 

It is interesting to consider further why the noise in 
the standard (mixed-order) Coulomb gauge calculation 
is problematic. Examining Eq. (24), we see that the 
backreaction (which depends on velocity v in this ex- 
pression) contains some perturbative amount of the "p- 
noise" , which is true thermal noise, obtained in our calcu- 
lation. Similarly, inspection also shows that the noise in 
the standard (mixed-order) Coulomb gauge calculation 
contains some perturbative amount of "p-backreaction" , 
which is the resistive damping that accompanies ther- 
mal noise, obtained in our calculation. In other words, 
the (non-thermal) noise in the standard (mixed-order) 
Coulomb gauge calculation contains some backreaction, 
and the backreaction contains some noise. Therefore, 
even in the classical limit of (24) we would have to enforce 
an artificial constraint upon the p-backreaction implying 
that the backreaction in the Coulomb gauge includes a 
non-zero average value of the p-noise. 

The physics may be clearer from another perspective. 
Recalling that rav = p — eA such that the velocity v 
depends on both the particle and field canonical coordi- 
nates, we see that velocity driving entails that the field is 
driven by both the system (particle) and itself. In other 
words, velocity driving of the field implies a perturbative 
feedback loop in the environment dynamics. Essentially, 
the field can self-generate stronger and stronger field ex- 
citations, albeit at an order beyond which the theory is 
accurate, and this can result in runways. In contrast, in 
our new equations of motion, Eq. (55), the field is driven 
by the canonical momentum p rather than the velocity 
v, and such pathological processes do not occur. 

Our next task is to demonstrate that the equations of 
motion are nonperturbatively stable. When combined, 
relations (57) and (58) reproduce the same damping 
as the Ford-O'Connell equation (48), which is already 
known to be stable, and the noises are perturbatively 
consistent. We can gain additional insight, however, by 
comparing our analysis to the analysis leading to (48). 
In this paper, we begin with an effective Hamiltonian, 
and work consistently to C(l/c 3 ) . With regard to the 
self-force this is essentially equivalent to taking both the 
dipole approximation and neglecting the A 2 term in the 
Hamiltonian. At C(l/c 3 ), we find that there is no con- 
straint on the cutoff and the backreaction is fully in- 
sensitive to particle structure. In contrast, in (48) the 
dipole approximation is made to the full Hamiltonian 
including the A 2 term, and instead a particle-structure 
form factor is assumed, which gives an effective cutoff in 
the particle-field interaction. The equations of motion 
are then order reduced, effectively to C(l/c 3 ), and the 
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resulting structure-independent equations of motion are 
found to be runaway free. The ultimate agreement be- 
tween these approaches is consistent with the fact that 
the contribution to backreaction from the A 2 term in the 
Hamiltonian is effectively negligible in the regime where 
the two methods of approximation are valid. 

Although it is already known that Eq. (48) is run- 
away free and causal, we now re-analyze these proper- 
ties within the present framework to see how this physics 
arises from with a purely effective theory framework. We 
will show that the dynamics of the open system are dis- 
sipative and stable in a manner analogous to the QBM 
analysis given in App. IV A. Consider a single particle 
and denote the system Hamiltonian 

P 2 

H sys = o Ve ^( x ) ■ 

An energy constraint can be obtained from either the 
Heisenberg equations of motion for H sys (£) or by inte- 
grating the (classical) Langevin equation (69) along with 
the second (70). Discarding the irrelevant transient terms 
(which can be shown to be bounded and runaway-free) 
we obtain the relation 

H sys (t) = H sys (0) + H 7 (i) + H e (t) , 

where 

H 7 (t) = -2^ 70 fdt'p(t'f (63) 
171 Jo 

is the energy lost to damping and 

(i) = - fdt'\ {i{t') ■ p(f ) + p(f ) • £(*')} (64) 

is the work done by the noise £(tf). The contribution 
from damping is manifestly a negative quantity. The 
noise is random and may do positive or negative work, 
but the damping only removes energy from the system 
(and delivers it to the environment and interaction). At 
least at this order, the damping is strictly local and so 
energy is lost to dissipation in a strictly uniform manner. 

It is important to note that the system here is given 
by the canonical variables (x, p), and H sys does not cor- 
respond to the mechanical energy of the particle, for the 
same reason that p 2 /2m is not the mechanical kinetic 
energy. From Eq. (20) the velocity and momentum differ 
by the vector potential, which includes both backreac- 
tion and noise. If the system momentum p relaxes under 
dissipative motion, then so does p and by extension the 
backreaction 70 p. Given noise, the system velocity fluc- 
tuates around the average 

(mv) € = (p + 2^-70 p^ , (65) 

implying that the system velocity is damped. If no ex- 
ternal forces are applied, the canonical and mechanical 
momenta approach each other (on average) in the late- 
time limit. 



B. Mass Renormalization 

In our Langevin equations (57)-(58), the p 2 /2m mass 
renormalization takes the form of magnetostatic self- 
interaction, which is ordinarily discarded in classical elec- 
trodynamics. Here we examine the counter terms in- 
volved in the renormalization and contrast them to the 
standard mass renormalization previously discussed in 
Sec. IIIC. 

Consider most simply the single particle theory in the 
dipole approximation. The resultant open-system equa- 
tions of motion are then given by 

i( t ) = ^W +2-^_ 70 p(t) (66) 

771 rcn ^barc ^barc 

p(t) = -eV0(x) , (67) 
in terms of the renormalized mass 

— = — +2^7(0). (68) 

'"Ten ^-barc ^barc 

Consistent with this order of perturbative analysis, we 
may express our Langevin equations as 

i( t )= + 2-^-70 p(t) --?-£(*), (69) 

772 rcn ^rcn ^rcn 

p(t) = -eV0(x) . (70) 

In this perspective, the renormalization of the p 2 /2m 
mass and the e/m (p • A) mass enter at different orders. 

It is well known that the instability of the Abraham- 
Lorentz equation arises from its implied negative bare 
mass for the system, which in-turn comes about if the 
high frequency cutoff A (or reciprocal radius c/r) ex- 
ceeds the characteristic frequency r" 1 [7, 14]. Yet to 
lowest order in r m (equivalently 1/c 3 ), the dynamics of 
charged-particle motion is insensitive to the high-energy 
details of the theory and is thus not problematic in that 
regime [7] . It is therefore not surprising that our pertur- 
bative approach yields cutoff-insensitive behavior, how- 
ever, the manner in which cutoff sensitivity is avoided is 
interesting. Whereas the standard radiation-reaction cal- 
culations involve a mass renormalization of mv 2 /2 given 
by Eq. (68), which runs the bare mass to negative in- 
finity in the high cutoff limit, our calculation runs the 
bare mass to positive zero in the high cutoff limit, and 
no pathological behavior is induced for any finite cutoff. 



C. Relativistic Kinematics 

For consistency in our quasirelativistic expansion of 
1/c, we should include the relativistic corrections to the 
single particle kinetic energy, as is standard in the Darwin 
Hamiltonian. Expanding the relativistic kinetic energy 
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gives 



H 



rel 



(mc) 2 + (p — e A) , 
, (p-eA) 2 (p-eA) 



2m 



8m 3 c 2 



+ ■ 



(71) 
(72) 



keeping terms of order C(l/c 2 ). In the absence of ex- 
ternal fields, this generates a second-order correction to 
Eq. (57) giving 



*(*) = 



1 - 



1 



P(*) 
mc 



P(*) 

m 



•2-j-ybp(t)--C(t), 
m 

(73) 

and otherwise the Langevin equations are identical. As 
in the Darwin Hamiltonian, this relativistic correction 
must be considered perturbatively. At the present order 
of perturbation theory, however, we can resum the term 
into the free velocity 



1 - 



1 



— V 

mc/ 



J3 

m 



O 



1 



\ mc I 



which is more convenient and better behaved in the equa- 
tions of motion. This can also be done in the effective 
Hamiltonian to ensure that the energy spectrum has a 
lower bound. 

Note that even nonperturbatively these kinematic 
corrections included in the stability analysis given in 
Sec. IV A will still yield dissipative backreaction, as all 
of these terms commute with the A • p interaction at 
the relevant order. The only modification will be in the 
definition of the canonical system Hamiltonian (now rel- 
ativistic), and the noise average of the velocity will be 
given by 



(m v) 



—Y 

mc / 



2— 7o P 

771 , 



(74) 



D. Electromagnetic Damping is Relativistic 

The standard Abraham-Lorentz equation (1) is com- 
monly referred to as "nonrelativistic" . Within the frame- 
work of a 1/c expansion, this equation is more accurately 
described as quasirelativistic. From this perspective, the 
damping force, which is 0(l/c 3 ), is intrinsically a rela- 
tivistic correction to the particle dynamics. Let us reex- 
press the (v/c) 2 magnetostatic and 1/c 3 damping forces 
in Eq. (57) as both arising from dynamical generators 



x, = 



d 
dt 




— P l 7o — P 

m, m. 



(75) 
(76) 



By comparison with Eq. (61)-(62), if the magnetostatic 
potential V is considered (v/c) 2 , then the damping gen- 
erator r is of relative order (r/c)(d/dt)(v/c) 2 . Thus the 
damping force can be interpreted as a more dynamical 
and more relativistic correction. Given that both of these 
generators arise from the same integration kernel with- 
out the full dipole approximation (see App. B), a higher- 
order "multipole" expansion in 1 j c can expect to include 
many more such terms of higher order. The relativis- 
tic nature of the expansion terms will be even clearer at 
higher orders. 



V. SUMMARY AND DISCUSSIONS 

We have derived new stochastic equations of motion 
(57)-(58) & (73) for multiple charged particles in the 
electromagnetic field. These equations of motion incor- 
porate the known relativistic corrections to the electro- 
dynamics of spinless charged particles to order C(l/c 3 ), 
including the electrostatic, magnetostatic, electromag- 
netic damping forces, and field fluctuations. More- 
over the equations of motion are manifestly causal and 
runaway- free. Our analysis shows that a 1/c expansion to 
the Coulomb-gauge Hamiltonian describes consistent and 
well-behaved nonequilibrium electrodynamics for spinless 
charged particles. Our consistent-order equations of mo- 
tion have a close correspondence with the order-reduced 
Ford-O'Connell equation, which is known to be well 
behaved. Whereas traditionally the Abraham-Lorentz 
equation has been considered an "exact nonrelativistic" 
equation, from the perspective of our analysis, patholo- 
gies in the standard Abraham-Lorentz equation are asso- 
ciated with inconsistent, mixed-order in 1/c, approxima- 
tions. Our view is that radiation reaction is intrinsically 
relativistic, though the Abraham-Lorentz equation may 
only fully capture its effect to lowest order in 1/c, and 
the order reduction used in deriving the Ford-O'Connell 
equation actually serves to rid the dynamics of inappro- 
priate mixed-order contributions. 

At C(l/c 3 ), mass renormalization is identified with 
the magnetostatic self-interaction and, at this order, the 
bare mass is positive for all cutoffs and the backreaction 
is fully insensitive to particle structure. An interesting 
question for future research is whether the cutoff insen- 
sitivity is preserved at higher-orders in the 1/c expan- 
sion, as to our knowledge no rigorously-derived relativis- 
tic equations of motion exist in the literature. For a sin- 
gle particle, we see that only at C(l/c 4 ) does the (dipole 
approximation to the) A(x) 2 interaction play a role, but 
at the same order multipole terms in the p • A(x) in- 
teraction must also be included for consistency. The 
standard results show that inclusion of the A(x) inter- 
action plus dipole approximation results in pathological 
equations of motion. It will be very interesting to see 
in detail how including the C(l/c 4 ) multipole terms in 
p ■ A(x) might resolve these pathologies and, in particu- 
lar, whether there will be runaway-free behavior for any 
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cutoff (i.e., full particle-structure insensitivity) or if a fi- 
nite (or bounded) cutoff will also be required. 

In conclusion, our results show that a 1/c effective- 
theory expansion provides useful new insights into 
charged-particle backreaction, and provides a systematic 
and consistent framework for extending to higher-order. 
This is a first step towards the goal of a consistent per- 
turbative approach to nonequilibrium relativistic electro- 
dynamics for charged particle motion, considered within 
a nonequilibrium yet Hamiltonian framework that in- 
corporates a well-defined description of stochastic noise. 
For instance, these results can be applied to analysis us- 
ing master-equation [ : 5] and influence-functional [16] for- 
malisms. Future work within this framework should also 
incorporate spin degrees of freedom. Relativistic effects 
like particle creation, however, will probably be more 
naturally described by describing the charged particles 
with the Dirac field. The advantage of this formalism, 
in contrast, is that it more naturally describes particle 
trajectories. 



The system + environment Hamiltonian is then given by 



f°° If 1 
H QB M = H sys + J dk - |(7r fc - eg k x.) 2 +u^q|j 



H sys = — + C/(x) , 



(A6) 
(A7) 



where, as determined by the gauge of our Lagrangian, 
p is the system momentum conjugate to x and 7r is the 
field "momentum" conjugate to q. 

Note that for m > and U (x) sufficiently well behaved, 
Hamiltonian (A6) is bounded from below in its energy 
spectrum. Therefore, under these conditions runaway 
solutions will not occur when the environment is initially 
described by a thermal state. 

Additionally note that the "bare" system potential in 
Eq. (A6) is given by 

C/ barc (x) = C/(x)+fe 2 / dfc^)x 2 , (A8) 
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and that the system + environment Hamiltonian can also 
be expressed as 



Hqbm — Hsys° — ex • 7T 



/•oo y 



in terms of the collective field operator 



7T = / dkg k -K k . 



Appendix A: Quantum Brownian Motion 

We begin our discussion with the Quantum Brownian 
Motion (QBM) Lagrangian which we have adapted in 
form and notation to better mirror the problem of backre- 
action in the electromagnetic field. This Lagrangian de- 
scribes a quantum system bilinearly coupled to a bosonic 
bath of harmonic oscillators and is traditionally used to 
model ordinary motional damping in quantum mechan- 
ics. 



£qBM — Ays + Ant + ^env j 
1 



-mx — U(x) , 



£ int =ex-Q, 



f°° 1 

£ c „v = J dfc-{q£-wgq!} , 



(Al) 
(A2) 
(A3) 
(A4) 



where x denotes the system position, <\ k denote the field- 
mode "positions" , and Q is the collective field operator 



The resulting Heisenberg equations of motion then dic- 
tate that the system is driven by the field 

x = - , (A9) 
m 

p = -W bare (x) + e7r, (A10) 
whereas the field modes are driven by the system 



qfe = TTfc - eg/cX , 
7T fc = -u\ q fe . 



(All) 
(A12) 



Solving for the field-mode evolution as driven by the 
system, we obtain the homogeneous + driven solution 

*k (t) = £ k (t) + e g k wl (G k *x) (t) , ( A13) 

£ fc (t) = G k (t) ir k (0) + G k (t) q k (0) , (A14) 

(A15) 



where the * product denotes the Laplace convolution 



Q = / dkg k c[ k . 



(A5) 



(A*B)(t)= dt'A(t-t')B(t'). (A16) 
Jo 
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The time-evolving field operator is then given by 



The damping kernel may then be expressed 



7T(t) = g(f)- 2e(/i*x)(t) , 

noise dissipation 

dkg k £ fe (i) , 



/*(*) = - / 



dk%±G k (t) 



(A17) 

(A18) 
(A19) 



where fi(t,t') = (j,(t — t') is the stationary dissipation 
kernel and £(f) is a Gaussian stochastic process for the 
initial conditions we assume: a factorized state of the sys- 
tem and environment, with the environment in a thermal 
state. 

Next we introduce the related damping kernel 



d 



7(M') 



<9i' 

» 2 

dfc^cosK(<-t')] 



(A20) 
(A21) 



which is necessarily positive definite and independent of 
the (factorized) initial state of the environment. The 
backreaction can then be expressed 

(Ai*x)(i) = ( 7 *x)(t)+ 7 (t,0)x(0) - 7 (t,t)x(t) , 

damping 



dissipation 



slip 



rcnormalization 

(A22) 

in terms of the positive-definite damping and where we 
have labeled the terms corresponding to the renormaliza- 
tion and initial short-time slip dynamics associated with 
the factorization of the initial state (see Sec. A3). 

Substituting our field solutions into the system equa- 
tions of motion, we obtain the quantum Langevin equa- 
tion 

mx(f) + 2 e 2 (7 * x)(i) + VC7(x) = e £(f) - e 2 7(f) x(0) , 
which reduces to 

mx(i) + 2 e 2 (7 * x) (t) + Vt/(x) = e £(i) , (A23) 
after the transient slip is taken into account. 

1. Ohmic Coupling and Local Damping 

Considering the damping kernel, which is given by 

7 (t) = / dk 9 ^ cos^f) . (A24) 
Jo 1 

If assume U)k = ck and gk ~ g up to some high-frequency 
cutoff A, then we may evaluate the integral as 

/ s g 2 f A , / \ g 2 sin(Ai) 
7 (f) = du cos(ujt) = - — - . (A25) 



2c 



7 (f) =2 7o 5 A (t). 

„2 



7o 



4c ' 

f A (*) - ^ 
Ttt 



(A26) 
(A27) 

(A28) 



in terms of the Dirac delta (t). In the high-frequency 
limit, the damping contribution to the Langevin equation 
becomes 



lim (7*x)(f) = 7ox(i), 

A— >oo 



(A29) 



or local damping. 



Renormalization 



For Ohmic coupling or local damping the quantum 
Langevin equation described by Eq. (A23) is phenomeno- 
logical, in the sense that its various parameters corre- 
spond to the physically measurable parameters at low 
energy. Assuming the Langevin equation to be phe- 
nomenological, note the bare system potential in Hamil- 
tonian perspective (A8) & (A10) as compared to the phe- 
nomcnological system potential E/(x) is 

C/ b arc(x) = [/(x) + 2e 2 7 (0)x 2 , 

2 

where 7(0) = |-A for local damping with a hard cutoff 
regulator. The renormalization is a quadratic term, re- 
gardless of whether or not the original model contained 
such a term. The QBM model typically proceeds from 
an x • Q interaction, where this renormalization does not 
naturally result from the Lagrangian theory. 



3. Factorized Initial Conditions 

If the operator noise £(i) in our Langevin equation is 
to be sampled from a thermal distribution which is (ini- 
tially) statistically independent from the system, then 
the initial state of the system and environment must be 



sys+cnv 



2c t 



a product state of the form p, 

a product of marginal phase-space distributions in the 
classical regime, and with the environment initially in 
equilibrium. This is an important simplification in our 
(and most other, e.g., [17-19]) analyses of the noncqui- 
librium dynamics of open systems. 

The consequence of assuming an initially uncorrelated 
system and environment must be carefully examined 
when studying radiation reaction, however, since acausal 
behaviors arise during the same very short time scale 
where the unphysical nature of a factorized state is rel- 
evant. It is therefore an important aspect of our anal- 
ysis that we are also able to apply recent results [10] 
showing that for classical or high-temperature electro- 
magnetic noise {hw sys <C k^T in Eq. (A32)) the initial 
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evolution of factorized states (or distributions) quickly 
leads to physical, dressed particle states without reintro- 
ducing the pathologies or instabilities in the dynamics 
that our analysis is intended to avoid. In the semiclas- 
sical or quantum regime, use of properly-correlated ini- 
tial states can mitigate the unphysical aspects of assum- 
ing initially factorized states entirely, without otherwise 
spoiling the results in this paper [20]. 



a. The Slip 

The transient slip in our Langevin equation is an 
initial-time pathology associated with vanishing corre- 
lations in the factorized initial conditions despite non- 
vanishing interaction strength between the system and 
field. In addition to the slip, there is a diffusive initial- 
time pathology, called jolts, which arise from correlation 
with the (quantum) zero-point fluctuations of the envi- 
ronment. The slip in particular was thoroughly analyzed 
in [10], where it was pointed out to generate the linear 
dynamical map 



i2e 7qx 



(A30) 



which maps states in a unitary fashion and preserves all 
kinematic moment invariants [21], including the uncer- 
tainty function. Therefore one can identify the post- 
slip state as a "renormalized" initial state which is more 
properly correlated with the environment and the prc- 
slip state as a "bare" initial state. If one only considers 
the classical regime, then jolting is not severe due to the 
lack of zero-point fluctuations in the environment. More- 
over, for a classical zero-temperature environment there 
is no noise causing any diffusion. Thus for this case one 
only needs to consider the renormalized initial states, ef- 
fectively discarding the slip term entirely. 



4. The Fluctuation-Dissipation Relation 



The Gaussian process has the noise kernel 

KMO = = ^-0, (A3i) 

which is stationary and positive definite for any station- 
ary initial state of the environment. For an equilibrium 
initial state of the environment the noise kernel is re- 
lated to damping kernel by the (quantum) fluctuation- 
dissipation relation (FDR) 



v{uj) = 7(0;) fiweoth 



2k B T J ' 



(A32) 



where j(u>) = f_ dt e~ ZUJt j(t). Essentially, the damp- 
ing kernel and temperature completely characterize 
Gaussian, thermal noise. 



Also note that as the coupling and environment are dy- 
namically linear, the damping kernel, being determined 
by the commutator, is independent of the state of the 
environment and it is the same whether in the classical 
or quantum regimes. In the classical regime we have the 
limit 



lim v(ui) 
?l->o 



%u)2k B T, 



(A33) 



and the classical fluctuations vanish in the zero temper- 
ature limit. In this limit (the classical vacuum) we can 
neglect the stochastic process £(t). 

In the quantum regime, the anti-commutator expecta- 
tion value (A31) is not sufficient to describe the statistics 
of the operator- valued stochastic process £(t). One ad- 
ditionally requires the commutator expectation value 



fi(t,t') 



1 

2th 



(K(tU(f)]> £ = M*-0 



(A34) 



which is given by the dissipation kernel, a state- 
independent quantity. Here the dissipation kernel is not 
generating backreaction, but consistent time evolution 
for the non-commuting stochastic process. The full quan- 
tum correlation is therefore given by 

a(t, f ) = (K(t) TL{t')) cm = {€(«) «f)> e , (A35) 
= v(t-t')-\-thfx{t-t'), (A36) 

where 7T(£) denotes the interaction-picture or Dirac- 
picture field operator and not the Heisenberg-picture field 
operator which we have already denoted 7T(t). 



5. Stability Analysis 

We will now show that the dynamics of the system are 
dissipative and stable under the very same conditions for 
which the system + environment Hamiltonian (A6) has 
a lower bound in its energy spectrum. Let us denote the 
canonical system Hamiltonian 



H 



sys 



2m 



t/(x) 



(A37) 



One may then calculate an energy constraint from either 
the Heisenberg equations of motion for H sys (£) or by in- 
tegrating the (classical) Langevin equation (A23) along 
with velocity. Accounting for the slip in our initial state, 
which only produces a finite change in energy, we obtain 
the relation 

H sys (t) = H sys (0) - H 7 (t) + H c (t) , (A38) 

in terms of the energy lost to damping H 7 (t) and energy 
generated by noise H^(t) 

H 7 (t) = e 2 fdt' fdt"j(t',t") l - {x(i') • x(t") + Mi!') ■ x(t')} , 

(A39) 



JO 
t 



H e (t) = efdt' - {£(*') • x(f ) + x(i') • m} ■ (A40) 
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The contribution from damping is a negative quantity 
as the damping kernel is a positive-definite kernel in a 
quadratic form. The noise is random and may drive the 
system erratically, but the damping may only remove 
energy from the system (and deliver it to the environ- 
ment and interaction). Therefore it is imperative that 
H sys have a lower bound in its energy spectrum. For our 
model, this necessarily implies that the system + envi- 
ronment Hamiltonian (A6) also has a lower bound in its 
energy spectrum. If this is the case then true runaway 
motion cannot occur. In the classical-vacuum limit, en- 
ergy is continually siphoned from H sys (i) until all motion 
ceases. 

Locally-damped energy is additionally simplified to 

H 7 (i) = 2e 2 7o x(i) 2 , (A41) 

which monotonically dissipates energy in time. Nonlo- 
cal damping can produce an instantaneous increase in 
system energy, though the cumulative effect is always 
dissipative. 



The driven solutions of each field mode are therefore 
given by 



a k .««(f) = e-—'a k .«,(0) (B2) 



In calculating A(x,) (10) we require evaluation of the 
field modes at the system location 



Appendix B: Explicit Calculation of Driven 
Quantum Field 

Most simply let us consider the field degrees of freedom 
as driven by one particle, using Hamiltonian (50). 



3 + ik - x ^a k , £fc (t) = \ {e+* k ' x ^,a k , £fc (t)} , (B3) 



aw, 



-iLV k ak, £fc 



E2m~ r _ 
, \ e 



' k x ,(e*-Pi)} ■ 
(Bl) 



which we have placed into symmetric form by commuta- 
tivity of the system and field operators. We now must 
consider the driven mode 



e +t kx ' (t) a k , 6fc (t) = \ ( e - 1 ^ t-k-x- (*)) , ^(q) } _ % y f dt' 1 e -«<-**-*-** (*)) , f e -H(o,*f -k^ (f )) ) £fe . p . ^ U 

(B4) 



which we have also placed into a manifestly Hermitian 
form. Note that any function /[k • x(i')] commutes with 
e k'P(t') given the Coulomb gauge constraint k-e^.. There- 
fore we may move €f. ■ p(i') to the most suitable side 
of e +lk x ( t For the velocity source, one can first de- 
compose the velocity into momentum • p(t') and field 
Ek ■ A[x(£')] and then note that by the previous argument 
both terms commute with any function /[k • x(i')]. For 
the position source, commutativity follows trivially. 



I 

In any case, one can produce relations such as (55) by 
substitution of the above expression into Eq. (10): 

a(x^) = t?{t)-y |l {(^ * PJ )(<) + (^ * Pj y(t)} , 

3 3 

(B5) 

where the dissipation convolution is given by 

04*p,-)(t)= / dt'^(t,t') P] (t'), (B6) 
Jo 

and with the noise and dissipation kernel microscopically 
determined to be 
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e fc 1 



\/2e cj k 2 



({-' 



(w fc t-k-x,(t)) 



, a k ,, 



fc (0)} + {< 



+i(w k t-k-x i (i)) t 



,(0,}) , 



(B7) 



1 / ^k^T EfcC ^ ? /p-*("fc*- k ^i(*)) p+tCwfci'-k.x^t')) _ +J(w fc t-k- x< (t)) -tCwfci'-k-x^i'))! mo^ 

(2tt) 3 / 2 J 2^ 2e Q uj k 2 I J ' [ 8J 



T 

2e w fc 



sin 



[w fc (t-0 - k • ( Xi (t)- Xj (t')) + 0(fi 2 fc 3 )] . 



(B9) 



In the classical calculation the phase factors in Eq. (B8) 
commute, however in the quantum-relativistic regime 
they give rise nontrivial interference terms as compared 
to the approximate (B9). The O(k) phases in (B9) con- 
tribute to the 1/c 3 dissipative forces (given microscopic 
structure) and 1/c 2 nondissipative forces (given multi- 
ple particles). These phases are discarded in the dipole 
approximation. The next order 0(h 2 k 3 ) quantum phase 
corrections are determined by two-time commutators of 
the system trajectories and are negligible to the order we 
work at. 

By 0(h) we only mean to keep track of the fact that 
these corrections contain commutator dependence, and 



vanish for the classical equations of motion. Planck's 
constant is not an expansion parameter which we con- 
sider in this work. The effect of the 0(k n ) or 0{l/c n ) 
phase corrections admit a simple dimensional analysis, 
as only the full-time trajectories x(t') are input into 
these highly-oscillatory integrals. Dimensionless correc- 
tions can therefore only be formed by derivatives of x(t) 
and factors of 1/c, though not necessarily v/c. Essen- 
tially the 1/c expansion appears to be a kind of multipole 
expansion here, which generalizes the usual dipole limit. 

From Eq. (B8) the positive-definite and Hermitian 
damping kernel is then given by 



_J f rf3fr V 6feC fe 1 / p -t(uJ*t-k-x*(t)) -H(u) fc t'-k-x 3 -(t')) , -H(w fc t-k.xi(t)) -i(w fc t'-k-xj (t'))\ /r>i n s 

?tt)3/ 2 J ^ 2e oW 2 2 I +e J ' [aW) 



T 



cos 



[cjfc(t-t') - k • (xi^-XjCO) + o(hk 2 )] 



(Bll) 



We finally note the relation 
d 



^■(M')=-|77^(M')+o(y) , 



(B12) 



and so, to the order we consider, we may integrate by 
parts the dissipation kernel into the damping kernel. The 
relativistic correction here does not precisely involve the 
dissipation kernel, but a kernel of with the same param- 



eter scaling and frequency sensitivity. 

The operator noise processes are also more complicated 
quantum mechanically, as even for a Gaussian state of 
the environment the noise cumulants cannot be evalu- 
ated exactly. However, to lowest order in 1/c, with the 
0(k) phases treated as quasistationary, the noise pro- 
cesses are dual to the 1/c 3 damping, with which they 
satisfy a fluctuation-dissipation relation. 



[1] J. D. Jackson, Classical Electrodynamics, 3rd ed. (Wiley, 
1998). 

[2] A. D. Yaghjian, Relativistic dynamics of a charged sphere: 
Updating the Lorentz- Abraham model (Springer, Berlin, 
1992). 

[3] F. Rohrlich, Am. J. Phys. 65, 1051 (1997). 
[4] R. Medina, J. Phys. A 39, 3801 (2006). 
[5] P. M. V. B. Barone and A. O. Caldeira, Phys. Rev. A 
43, 57 (1991). 



[6] G. W. Ford, J. T. Lewis, and R. F. O'Connell, Phys. 

Rev. A 37, 4419 (1988). 
[7] G. Ford and R. O'Connell, Phys. Lett. A 157, 217 

(1991). 

[8] Dalibard, J., Dupont-Roc, J., and Cohen- Tannoudji, C, 

J. Phys. France 43, 1617 (1982). 
[9] B. L. Hu and A. Matacz, Phys. Rev. D 49, 6612 (1994). 
[10] C. H. Fleming, A. Roura, and B. L. Hu, Ann. Phys. 326, 
1207 (2011). 



17 



[11] J. T. Hsiang, T. H. Wu, and D. S. Lee, Foun. Phys. 41, 
77 (2011). 

[12] G. Breit, Phys. Rev. 36, 383 (1930). 
[13] G. Breit, Phys. Rev. 34, 553 (1929). 
[14] E. J. Moniz and D. H. Sharp, Phys. Rev. D 15, 2850 
(1977). 

[15] C. H. Fleming, N. I. Cummings, C. Anastopoulos, 
and B. L. Hu, "Non-Markovian entanglement of two- 
level atoms in an electromagnetic field," (2010), 
arXiv:1012.5067 [quant-ph]. 

[16] C. Anastopoulos and A. Zoupas, Phys. Rev. D 58, 105006 
(1998). 



[17] R. P. Feynman and F. L. Vernon, Ann. Phys. 24, 118 
(1963). 

[18] A. O. Caldeira and A. J. Leggett, Physica A 121, 587 
(1983). 

[19] B. L. Hu, J. P. Paz, and Y. Zhang, Phys. Rev. D 45, 
2843 (1992). 

[20] C. H. Fleming, A. Roura, and B. L. Hu, "Initial 
state preparation with dynamically generated system- 
environment correlations," (2011), arXiv:1101.2668 
[quant-ph] . 

[21] A. J. Dragt, F. Neri, and G. Rangarajan, Phys. Rev. A 
45, 2572 (1992). 



